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Abstract. The main production mechanism for supersymmetric Higgs particles at hadron colliders 
crucially depends on the size of their Yukawa couplings to bottom quarks. For sufficiently large 
tan P the total cross section for some of the neutral Higgs bosons in the MSSM is dominated by 
bottom-quark fusion. After an introduction to bottom-associated Higgs production, we discuss the 
complete 0{a) electroweak and 0{as) strong corrections for the bb-fusion channel in the MSSM. 
Choosing proper renormalization and input-parameter schemes, an improved Born approximation, 
constructed from previously known results, can absorb the bulk of the corrections so that the 
remaining non-universal corrections are typically of the order of a few per cent. Numerical results 
are discussed for the SPS benchmark scenarios. 

PACS. 12.60.Jv supersymmetric models - 14.80.Cp non-standard-model Higgs bosons 



1 Introduction 

The Higgs mechanism is a cornerstone of the Standard 
Model (SM) and its supersymmetric extensions. Thus, 
Higgs bosons are intensively searched for at the up- 
graded proton-antiproton collider Tevatron, followed 
in the near future by the proton-proton collider LHC. 
In this talk, we concentrate on the precise prediction 
of the total Higgs-boson production cross section at 
the LHC based on the results in Ref. [T]. 

In the SM, the total production cross section for 
Higgs bosons H at the LHC is dominated by gluon 
fusion. Higgs radiation off bottom quarks [5] 



pp/pp-^bb0"-fA: 



(1) 



with (fp = H, is a negligible contribution. The relevant 
bottom Yukawa coupling A^^ is known to be small be- 
cause it is determined by the ratio of the small bottom- 
quark mass mb and the known vacuum expectation 
value (VEV) v of the SM Higgs field, A^^ = rrn,/v. 

In contrast, in the MSSM, bottom-associated pro- 
duction of neutral Higgs bosons, (/)° = h°, H°, A°, can 
dominate the total cross section at large tan/?. Two 
different Higgs doublets are needed to generate masses 
for up- and down- type fermions. These two Higgs dou- 
blets Hu and Hd acquire VEVs Vu and Vd, respectively, 
and one defines tan/3 = Vu/vd- While v'^ = + v'^ is 
fixed by the gauge-boson masses, the ratio tan/3 is a 
free parameter. For large tan/? the down-type VEV 
Vd is small and the Yukawa coupling of the down-type 
Higgs doublet is enhanced with respect to its SM value. 



The mass mb is not small due to a small Yukawa cou- 
pling, on the contrary, the relevant VEV Ud is small. 
For tan/3 ~ 0(50), the bottom Yukawa coupling is as 
big as the top Yukawa coupling in the SM. 

The couplings of the CP-even Higgs-boson mass 
eigenstates are determined by the mixing of up- and 
down-type Higgs fields characterized by the mixing an- 
gle a. For the Yukawa couplings to b quarks one finds 
in the MSSM 
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For sizeable mixing in the Higgs sector, the cross sec- 
tions cr for the b-associated production of all neutral 
MSSM Higgs bosons are enhanced for large tan (3, i.e. 
a cx tav? (3. However, one finds sin a — > — cos/3 if the 
mass M\ of the pseudoscalar Higgs boson, the sec- 
ond input parameter of the MSSM Higgs sector, is 
large compared to the Z-boson mass. In this limit, the 
lighter CP-even Higgs boson h" is SM like and shows 
no enhanced bottom Yukawa coupling. Nevertheless, 
the total cross section for the two heavy neutral Higgs 
bosons H" and A*^ is dominated by b-associated pro- 
duction. 

Current searches for bottom-Higgs associated pro- 
duction in the MSSM at the Fermilab Tevatron ex- 
clude values tan/3 > 50 for light « 100 GeV [3]. 
For a recent sensitivity study at the LHC see Ref. ^ . 

The theoretical description for bb (jp production 
can start from different initial states for the hard scat- 
tering process. The b quarks are either generated from 
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gluon splittings within the hard process or they are 
considered to be part of the proton, i.e. the gluon split- 
ting is factorized from the hard process. 

In the so-called four-flavour number scheme (4FNS) 
with no b quarks in the initial state, the lowest-order 
QCD production processes are gluon-gluon fusion and 
quark-antiquark annihilation, gg bb^° and qq — > 
bb (jP , respectively. In this framework the splitting of 
gluons into bb pairs is treated retaining the full depen- 
dence on the bottom mass. The complete kinematics 
of the 2-^3 process is available so that the bottom 
jets in the final state can be used for tagging and back- 
ground suppression. 

However, for a hard process involving a large scale, 
e.g. the Higgs-boson mass, the b quark is effectively al- 
most massless. Higgs production is thus dominated by 
events where gluons split into nearly coUinear bb pairs. 
Consequently, the inclusive cross section for gg — > bb (fP 
contains large logarithms hi{^p/mh), where the large 
scale fj,F ~ M^o corresponds to the upper limit of 
the coUinear region up to which factorization is valid. 
Hence, the perturbative expansion in lii{fip /mh)a.s will 
eventually break down for large Higgs masses and the 
perturbation series has to be reorganized. The loga- 
rithms hi{fip /myy) can be summed to all orders in per- 
turbation theory by introducing bottom parton den- 
sities. This defines the so-called five-flavour number 
scheme (5FNS) [5]. In this scheme, the leading-order 
(LO) process for the inclusive bb 0° cross section is bb 
fusion, 

bb ^ 0° . (3) 

The next-to- leading order (NLO) cross section in the 
5FNS includes 0{as) corrections to bb (jP and tree- 
level processes like gb h<tP . For developments on 
corrections to the latter subprocess see Ref. [6] . 

The use of bottom distribution functions is based 
on the coUinear approximation, i.e. outgoing b quarks 
are considered to have small transverse momentum 
and to be part of the proton remnant. There is no 
theoretical control over additional b jets at LO. 

To all orders in perturbation theory the four- and 
five-fiavour schemes are identical, but the way of or- 
dering the perturbative expansion is different, and the 
results do not match exactly at finite order. However, 
numerical comparisons between calculations of inclu- 
sive Higgs production in the two schemes 7,8,..9lfT0] 
show that the two approaches agree within their re- 
spective uncertainties, once higher-order QCD correc- 
tions are taken into account. 

There has been considerable progress in improv- 
ing the cross-section predictions for inclusive associ- 
ated hh (jP production by calculating NLO-QCD [Tj 
[9] and SUSY-QCD TTj corrections in the four-flavour 
scheme, and NNLO QCD corrections |:12^ in the five- 
flavour scheme. In the 5FNS, the QCD scale uncer- 
tainties have been reduced to the 10% level such that 
radiative effects from the electroweak sector become of 
interest. 

The complete one-loop QCD and electroweak cor- 
rections for the decay of MSSM Higgs bosons to bot- 
tom quarks have been presented in Ref. |13| . Recently, 



complete supersymmetric QCD and electroweak cor- 
rections to the hadronic production cross section have 
been presented in Ref. [1]. These results and their re- 
lation to known universal corrections are discussed in 
the following sections. For technical details and deriva- 
tions we refer the reader to Ref. [1]. 

2 Radiative corrections 

In b-quark fusion, bb — s- (jP , there are universal ra- 
diative corrections which lead to the definition of the 
improved Born approximation for the partonic cross 
section 

sin^ Qcff ( 1 - ^b/ (tan (3 tan acs) \ ^ 

aiBA = asM ( l + ^b tanaeff/tan/g y 

cos^ P \ 1 + /Jb J 

(4) 

where 

o-SM = 5(1 - M^o/s) . (5) 

We denote the partonic CMS energy by Vs , M^o is 
the mass of the produced Higgs boson, is the Fermi 
constant, and rnh{fJ.B.) is the running bottom mass at 
the renormalization scale /Lt/j. Potentially large radia- 
tive corrections are encoded in the parameters Z\b and 
ttoff to be briefly explained in the following. 

While b quarks do not couple to the Higgs field 
Hu at tree level, this interaction is radiatively induced 
at the one-loop level, e.g. via the sbottom coupling 
to Hu in SUSY-QCD loops. This induces a shift 
in the relation between the b mass and the respective 
Yukawa coupling. This shift is proportional to tan/3 
and, thus, for large tan/3 the corresponding correction 
is sizeable. It has been shown [14] that the correction 
can be resummed and that it affects the cross section 
according to eq. ([4]). 

Radiative corrections can also have a large impact 
on the mixing of the Higgs fields to form the CP-even 
mass eigenstates h° and H". The bulk of these cor- 
rections can be absorbed in a loop-corrected, process- 
independent effective mixing angle acff which replaces 
its tree-level counterpart in eq. (jl]). 

Precise definitions of Z\b as well as acfi are given in 
Ref. [T], where we also describe in detail the calcula- 
tion of the complete SUSY-QCD and electroweak cor- 
rections including the renormalization of the MSSM 
Higgs sector. The bottom mass has been renormal- 
ized in such a way that the corrections due to Z\b 
are absorbed into the input value for mb. Thus, the 
numerical value for this effective mb quantifies the 
tan/3 enhanced corrections. This procedure automat- 
ically avoids double counting for the Ay, corrections. 
When we relate the corrections from the full calcula- 
tion to (TiBA we also carefully avoid double counting 
with respect to corrections already contained in acff. 
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MSSM Higgs-bosons from bottom quarks 
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Table 1. The eflfective bottom mass and the NLO MSSM cross section pp (bb) h^/H^/A^ + X at the LHC (^s = 
14 TeV) in the SPS 4 scenario including the cumulative corrections due to the different classes of corrections. See text 
for details on the different contributions. (Table taken from Ref. |Tj) 



3 Results 

All the results in this section are calculated in the DR 
scheme for tan/3. The renormalization and factoriza- 
tion scales are set to /i/? = M^o and fip = M^a/A, re- 
spectively. We use the MRSTQED2004 PDF jj which 
also allows the inclusion of the photon-induced par- 
tonic channels at NLO. The input-parameter scheme 
and the numerical input are specified in Ref. [1]. To 
further improve our NLO results, we use two-loop im- 
proved Higgs self-energies provided by the program 
package FeynHiggs 16J. 

Within the MSSM, we first focus on the radiative 
corrections and total cross sections in the SPS 4 bench- 
mark scenario (tan (3 — 50) [17] which was designed to 
give large cross sections for heavy Higgs bosons. At the 
end of this section we also show results for the other 
SPS points. While most of the SPS scenarios are in 
conflict with experimental data by now, they are still 
valuable because they cover typical SUSY scenarios in 
different regions of parameter space. 

In Table [U we show the effect of the various higher- 
order corrections on the effective b-mass. Starting from 
the running QCD mass at the scale of the Higgs-boson 
mass, the shifts from SUSY-QCD {A^) and the elec- 
troweak sector [A^'^'^^) are included. The correspond- 
ing cross sections a are first shown at NLO QCD. As 
can be seen, the QED corrections are generally very 
small after mass factorization. The summation of the 
tan /3-enhanced MSSM-QCD and MSSM-weak correc- 
tions has a significant effect on the cross sections for H" 
and A" production. The light Higgs boson h° is SM- 
like and the summation of terms oc tan (3 has thus no 
sizeable impact. Employing a loop- improved effective 
mixing angle is numerically relevant only for h" 
production. The cross sections in the last-but-one row 
of Table [T] correspond to the improved Born approxi- 
mation (TiBA dressed with QCD and QED corrections. 
The full MSSM cross sections, including all summa- 
tions and the remaining non-universal 0{as) and 0{a) 
corrections, are displayed in the last row of the table. 

The bulk of the MSSM-QCD and -weak corrections 
can indeed be absorbed into the above definition of 
CTiBA- The remaining non-universal corrections in the 
complete MSSM calculation turn out to be quite small, 
below approximately 2%. 

In Fig. [T] we show the impact of the complete super- 
symmetric 0{as) and 0{a) corrections defined relative 



to the improved Born approximation ctiba for different 
values of the on-shell mass M^. All other MSSM pa- 
rameters are kept fixed at their SPS 4 values. The size 
of the non-universal corrections does not exceed 3% 
for H°/A° production except for special model param- 
eters where the Higgs masses are close to the produc- 
tion threshold for pairs of sparticles. These unphysical 
singularities can be removed by taking into account 
the finite widths of the unstable sparticles. The size 
of (5mssm for h'^ depends very sensitively on the def- 
inition of the effective mixing angle employed in 
c^iBA- Note that in any case h" is SM-like at large Af^ 
so that h'' production in bottom fusion is most likely 
of no phenomenological relevance. 

It is important to emphasize that the non- universal 
MSSM corrections Jmssm at large tan (3 are quite sen- 
sitive to the choice of the b-mass input value within the 
one-loop corrections which is not fixed by the renor- 
malization procedure. There are terms that grow as 
TO^ tan^ /? which are not included in Z\b- For the SPS 4 
scenario the sensitivity on the b-mass input is shown 
in Fig. [21 The absolute size of the non-universal cor- 
rections varies between approximately zero and —6% 
for the phenomenologically relevant H°/A" produc- 
tion, depending on whether a massless approximation, 
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Fig. 1. Full MSSM corrections (5mssm defined relative to 
(JiBA as a function of the Af^ pole mass. All other MSSM 
parameters are fixed to their SPS 4 values. (Figure taken 
from Ref. [I]) 
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Fig. 2. Full MSSM corrections (5mssm defined relative to 
(Tib A as a function of the rrib input. The corrections for 
H" and a" lie on top of each other. (Figure taken from 
Ref. [1]) 



the effective running mass, or the pole mass is chosen 
as b-mass input. Although we assume that the running 
mass, including the corrections from Z\b (as used for all 
the shown results), is a sensible choice, the sensitivity 
of the NLO correction to the b-mass input constitutes 
a theoretical uncertainty which cannot be resolved at 
the NLO level. 

Table [2] displays the cross sections along with the 
non-universal corrections from the full calculation for 
the different SPS points. It shows that these residual 
corrections are small and do not exceed 2% for H°/A'^ 
production in a wide range of MSSM parameters. 



4 Conclusions 

We have given a brief review on Higgs production in 
association with bottom quarks focussing on the ba- 
sic concepts and precise predictions for Higgs produc- 
tion in bottom-quark fusion. The leading supersym- 
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metric higher-order corrections can be taken into ac- 
count by an appropriate definition of the couplings in 
an improved Born approximation. The remaining non- 
universal corrections are small, typically of the order of 
a few per cent. The theoretical uncertainty connected 
to the input value of the b-quark mass within the one- 
loop correction is emphasized. 
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